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1. Introduction 

Many problems in physics require to study how the overall properties of 
a system are influenced by its boundary. An archetypal task is looking for 
dependence of the lowest eigenvalue of the (minus) Laplacian, defined in a finite 
domain, on the shape of this domain. This problem is related to topics like 
the calculation of the ground state energy of a free quantum particle in a finite 
domain, the lowest cutoff frequency of the dominant mode in a waveguide, or the 
lowest tone of a drum (HQ. We shall formulate this mathematical problem in 
the context of the diffusion (probabilistic) theory, which is exploited often in the 
chemical physics: The geometry of a molecular boundary influences the reaction 
rates Q, or the dynamics of particles diffusing in confined systems 0, @. Our 
attention is focused on the survival probability of a particle diffusing in a finite 
domain, which is defined later. We show that this quantity interrelates the 
boundary characteristics and the lowest eigenvalue of the Dirichlet Laplacian in 
an interesting way. 

To introduce the notation, we consider a finite domain f2 of points r with 
a smooth boundary dQ. The spectrum of the (minus) Laplacian, say with the 
Dirichlet boundary condition (BC), is given by 

-A^(r) = A0(r) r€(], 

0(r) = real!. ( ' 

The eigenvalues < Ai < A 2 < A3 ••• < Xj < ■ ■ ■ form a discrete set@]. The 
corresponding eigenfunctions <f>\, (j>2, 03, . . . , <f>j, ■ . ■ form an orthonormal basis of 
real functions, 

/ dr^(r)^(r)=^; (2) 
the eigenfunctions satisfy the completeness relation 

5> J -(r)&(r , )=*(r-r'). (3) 
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In the diffusion theory, the conditional probability p(r,i|r o ,0) of finding a 
particle at a point r £ ft at time t > 0, if it started from ro £ f2 at to = 0, is 
governed by the diffusion (heat) equation 

d -^p^ = A P (r, t \r ,0), (4) 

where the diffusion constant is set to 1. This equation has to be supplemented 
by the initial condition p(r,t — 0|ro,0) = S(r — ro) and by the Dirichlet BC 
p(r, t|ro,0) = for r £ dil, which reflects absorption (disappearance) of the 
particle hitting the boundary. The conditional probability can be expressed in 
terms of the eigenvalues and eigenfunctions of the Dirichlet Laplacian as follows 

p(r,t|r ,0) =£>i(i-o)4>i(r)e- A '- t . ( 5 ) 

3 

The "local" survival probability is defined as 

S(t,r ) = f dr/j(r,t|r ,0). (6) 

It represents the probability that a particle, localized at a point ro £ O at 
the initial time to = 0, remains still diffusing in the domain 51 at time t > 
unabsorbed by the boundary. Since the particle inserted at the boundary at 
to = is immediately absorbed, S(t, ro) satisfies also the Dirichlet BC 

S(t,r o ) = forr £<9ft. (7) 

The quantity of our interest is the "global" survival probability S(t), defined as 
the average of the local one over the whole domain, 



(8) 



It represents the probability of finding the particle in Q at time t > 0, if it was 
distributed uniformly with the density l/|f2| over the whole domain at t = 0. 

The same analysis can be applied to the heat contained in the domain f2, 
heated at some nonzero temperature, whose boundary dQ is cooled to the zero 
temperature for t > 0. S(t) is then proportional to the heat content in the 
domain and this is the name under which this quantity is known in the mathe- 
matical literature d i, [l(J [TT], 0, El • 



First let us review some basic properties of S(t). At t — 0, it starts from 
iS"(0) = 1 as the particle just inserted somewhere in the domain is surely not yet 
absorbed by the boundary; the consequent equality Y]^ = 1 follows trivially 
from the completeness relation ([3]). The monotonous decay of S(t) at small t 
is determined by the geometry of the boundary, as only the particle inserted 
close to the boundary have a chance to be absorbed very soon. Like for any 
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spectral function, the small-t expansion of S(t) can be obtain explicitly. It is an 
expansion in \fi of the form 



oo 

S(t) = l + ]>>^/ 2 , (9) 

3=1 

for any dimensionality of the domain [l2l. Il3j|. The coefficients {<7j} were shown 
to satisfy a recurrence scheme, containing certain differential operators acting on 
the shortest distance l(r) of points r from the boundary 917 [lj; EH- The forms 
of the coefficients become soon very complicated, only first few of them are 
available explicitly. The first aim of this paper is to approximate reasonably the 
coefficients {<Xj} in higher orders by local curvatures of the domain boundary. 

It is clear that the whole eigenspectrum contributes to the coefficients {<Jj}. 
On the other hand, according to ([8]), the long-time decay of the survival prob- 
ability to is governed by the lowest eigenvalue of the Dirichlet Laplacian Ai, 

S{t) ~ ~fl exp(-Ait) for t -> oo. (10) 

Thus the survival probability interconnects the geometric characteristics of the 
domain boundary, contained in the small-t expansion of S(t), with the overall 
characteristics of the system, represented by the lowest eigenvalue Ai in the long- 
time limit. Our second aim is to apply some appropriate interpolation scheme 
between the small-t and large-t expansions. In this way we can investigate how 
the lowest eigenvalue of the Dirichlet Laplacian in a domain depends on the 
geometry of its boundary. 

We are inspired by the method of t-expansions for studying the ground state 
energy of many-body systems 14]. The method is based on an energy functional 
E(t), depending on time t. This functional has an analytic small-t expansion in 
terms of connected moments {//}, E(t) — Yl'jLo(~ tYlj+i/j" The ground state 
energy Ai follows from the long time behavior Ai = lim t _ i . 00 E(t). There exist 
several interpolation schemes between the small-t and large-t expansions (lH [lil 



121 ll8( which enable one to calculate Ai in terms of connected moments {Ij}. 
The interpolation schemes are based on various combinations of exponential 
decays of E(t) to the asymptotic value Ai; the fact that the small-t expansion is 
analytic plays a fundamental role. The method of t-expansions was successfully 
used for calculation of the ground state energy of molecul es flil l2(i Ell. [22^ . or for 
the study of phase transitions in infinite lattice systems [18(. Since the small-t 
expansion of the survival probability @ is in powers of y/t, the interpolation 
schemes working in t-expansions are not straightforwardly applicable to our 
problem. 

As was indicated, this paper is devoted to and brings new results in two 
topics: 

• A reasonable approximation of the coefficients {crj} of the small-t expan- 
sion ([9]) in high orders. 

• The search for an interpolation formula between the small-t ([9]) and large-t 
(fT0|) expansions of S(t). 



3 



It turns out to be advantageous to perform the analysis in the Laplace space. 
In Sect. 2, we introduce the Laplace transform for the local and global survival 
probabilities and derive a differential equation for the former. In Sect. 3, we 
use the exact solution for a disk domain to approximate reasonably the coef- 
ficients {<Tj} of the small- i expansion © up to high orders which is desirable 
for the proposed interpolation method. The first four coefficients are identical 
to the exact formulas [lH EH, the next approximative ones involve correctly the 
powers of the boundary curvature integrated along dil. Sect. 4 is devoted to an 
interpolation method between the small-t and large-t expansions of S(t), based 
on a variant of the Pade approximant of the Laplace transform of S(t). Our 
algorithm is tested on the circular and elliptical domains. 



2. Laplace transform of S(t) 

Our method is formulated in the Laplace space. We introduce the Laplace 
transforms with respect to t, with s 2 as the Laplace parameter, for both the 
local survival probability 

/•OO 

r(s,r )=/ d£e- s2t S(*,r ), r(s, r ) = for r e dtt (11) 
Jo 

and the global survival probability 

r(s) = / dte~° 2t S(t) = £ 



(12) 



s 2 + Xj 



The Laplace transforms are related by 

r( S ) = p^dr r( S ,r ). (13) 

In the Laplace picture, the small-i behavior of S(t) ([9]) is represented as a large-s 
expansion of r(s) in 1/s: 

i r(j/2 + i) 

3=1 

The large-i behavior of S(t) is represented as a regular small-s expansion of t(s) 

9 

m s : 

oo oo 2 

T(,s) = 5>l) fe ^£^, (15) 

k=0 j=l 3 

with a dominant contribution of the lowest eigenvalues. 

The Laplace transform of the local survival probability satisfies a differen- 
tial equation which is derived by the following procedure[23j. If we apply the 
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Figure 1: Fermi coordinates for general domain Q with boundary dQ. 



conjugated Laplacian A + (acting upon the coordinates of ro) to both sides of 
(fTTj) and use the conjugated of the diffusion equation ([4]), we obtain 

pOO p p poo 

A+r( S ,r )= / dte" s2t / drA+p(r, t\r , 0) = / dr / dte^dtpir, t\r , 0). 
Jo Jn Jn Jo 

(16) 

Integration by parts in t then implies 

Ar(s,r) - s 2 t(s,v) = -1, (17) 

where we abandon the subscript and replace A + by A. This differential 
equation is supplemented by the Dirichlct BC 

t(s, r) = forredfi. (18) 

In the remaining part of the paper, we shall restrict ourselves to two-dimensional 
star-shaped domains with only one smooth boundary. To express r(s) for such 
domains, it is convenient to average the differential equation (|17[) over £1 and 
then to apply the Gauss-Ostrogradsky theorem: 

s 2 r(s) = 1 + THi / ^ Ar(.s, r) = 1 - ±- [ d J^^ . (19) 
l"| Jn m Jan °Z «=« 

Here, we introduced the Fermi coordinates (2^| (see Fig. 1); lo runs around the 
domain boundary and £ is orthogonal to lo such that £ = at the boundary dfl 
and £ > inside the domain fl. 



3. Approximation of S(t) by local curvature 

3.1. Disk domain 

The simplest star-shaped domain is a disk of radius R, with the boundary 
defined by r(ip) = R. The disk geometry is the crucial one in our approach. 
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The eigenfunctions of the Dirichlet Laplacian for the disk are expressed in polar 
coordinates as follows 



<Pn,i(r,ifi) = c ni iJi(ry/\ n j)cos(ltp), 



(20) 



where n = 1,2,... and I = 0,1,2,... are the radial and angular quantum 
numbers, respectively, c n j are the normalization constants and the eigenval- 
ues A rii ; = (zi, n /R) 2 are determined by the zeros z^ n of the Bessel functions 
of the first kind Ji{zi, n ) = 0. The lowest eigenvalue Ai = Ai.o = z^/R 2 = 
5.783196. . ./R 2 . Only the radial modes <f> n ,o(r,<p) — Jo(zo lTl r / R) j \/nRJi(zo, n ) 
give nonzero integrals J Q rdrd(p(j) n fi(r,(p) — 2^/ttR/ z$^ n in equation and 
therefore contribute to S(t), i.e. 



oo 



(21) 



This means that the survival probability does not carry information about the 
whole spectrum of the Dirichlet Laplacian, only the eigenstates with the angular 
quantum number I — are involved. This is a specific property of the disk 
domain. 

The representation (|2ip , written in terms of the zeros of the Bessel function 
Jo, is not suitable for our purposes. Instead, we use the differential equation 
(flTj) supplemented by the Dirichlet BC (fT8|). In polar coordinates, the general 
solution of (p~T|) . regular at r — > 0, reads 



r(s;r, tp) 



1 



sr)e 



Hip 



(22) 



Z=-c 



where {/;} are the modified Bessel functions of the first kind. Since r(s;r,ip) 
is real, the complex coefficients {77} must satisfy t_/ = r ; *. The boundary 
condition r(s; r = R, (p) = implies 



TO 



1 1 



Ti = for I ^ 0. 



Consequently, 



-{s;r,ip) = r(s,r) = \ 



1 - 



Io(sR) 



(23) 



(24) 



The averaging of this relation over the disk area leads to the explicit form of 
the Laplace transform of the survival probability [10| 



2/i (si?) 



sRI {sR) 



(25) 



To obtain the small-s and the large-s expansions of t(s), we use the corre- 
sponding expansions of the Bessel functions Iq(sR) and J1 (sR) [25| . The small-s 
expansion becomes 



In the limit s — > oo , using the asymptotic expansions of the Bessel functions Ig 
and I\, we obtain 



s 2 r{s) = 1 - — 
sit 



1 - 



2sR 



1 1 25 

(sR) 2 ~8(sR) 3 ~ 128(si?) 4 



(27) 



Comparing the series representations ([9]) and (|14p . we get the corresponding 
small-i expansion of the survival probability: 



S(t) 



1 



4v^ 



t 

R 2 



t 3 ' 2 



t 2 

8i? 4 



5t 5 / 2 
24V^i? 5 



(28) 



5.^. Convex domain 

Being motivated by the exact solution for the disk, we now propose an 
approximate form of the large-s expansion of t(s) for an arbitrary star-shaped 
domain. In the limit s — > oo, r(s,r) ~ 1/s 2 inside the domain and it decreases 
to 0, given by the Dirichlct BC, only in a thin layer near the boundary. This can 
be seen from the result ([24| if we use the asymptotic expansion Io{sr)/Io(sR) ~ 
p -s(R-r) f or s _> oo. Our strategy is to replace the adjacent part of the boundary 
adequately by some simpler curve, for which the function t(s; r) is known, and 
use this function instead of the exact one in the formula (|19p. 

The simplest way is to approximate the convex parts of the boundary by 
circles of radius R(uj), inverse to the curvature of the boundary at a given point 
u (see Fig. 1). Then, instead of the true t(s;w,£), we adopt the disk solution 
(|2"4"|) with r replaced by r = R(w) — £, so that 



<9t(s;w,£) 



_ h( S R{io)) 
£=o ~~ sI (sR(uj))' 



(29) 



Applying the asymptotic s — >■ oo expansions of the Bessel functions in this 
formula, the relation (1191) can be rewritten directly as an expansion in 1/s: 



s 2 f(s) 



1 - [ duj 
s\n\ Jan 

25 



1 



2sR{lu) 8s 2 R 2 {lu) 8s 3 i? 3 (w) 
13 



128s 4 i? 4 (a;) 32s 5 R 5 {lu) 



(30) 



Tilde on the top of f(s) means the approximation by the local curvature. Ac- 
cording to the representation (|14p . the corresponding coefficients {oj} are given 
by 



<?2 



OS 



1 1 f 



(31) 
(32) 
(33) 
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* - hwiL^' (34) 

" ^mL M <35) 
- - 55j5j /„**■«»>• < 36) 



etc., where fc(w) = 1/R(w) is the curvature at the point a) <E dCl. It follows 
from the construction that the coefficients can be expressed explicitly to very 
high orders. If the boundary is defined in polar coordinates, r = r((p) with 
(f € [0, 27r], the length element of the boundary dCl is given by 

duj = \Jr 2 + r ,2 dtp (37) 

and the curvature by 

*(") - r ;tf,-^" . (38) 

if —j— f'^ 

The integral in (|3~Tj) for is the length of the boundary |9f2| for an arbitrary 
domain, i.e. 

2 \dn\ , . 

81 - "7f W- (39) 

The substitution of the relations (157]) and (j3"8")l into (j3"2")l yields 

7r 1 / , r — rr 7r 



* 2 = W + 2Wn ^7M^ = M' (40) 



since the integral 

/■2tt 



dy- 







arctan ( — 
r 



(41) 



for an arbitrary r(ip). Note the universal form of |fi|(72 = tt valid for any star- 
shaped domain. The form of the higher coefficients depends on the particular 
shape of the domain boundary. 

3.3. Concave domain 

If the boundary possess some concave parts, these parts can be replaced 
by circles centered outside of the domain. Consequently, the function r(s;w,£) 
in (IT9]) is replaced by the solution of the differential equation ([T7]) for a thick 
annulus of the inner radius R — R(cu) and the outer radius going to infinity, 
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where Hu (z) = J v [z) + iY lJ {z) denotes the Hankel function^ 5]. After the 
derivation with respect to £ = r—R at the boundary £ = 0, the ratio iH^ (isR) / H^p [isR) 
can also be expressed using asymptotic s — ¥ oo expansions 25], with the result 



_ 1 



1 1 1 25 13 



2sR 8(sR) 2 8(sR) 3 128(si?) 4 32{sR) 5 

(43) 

In comparison with the formula (|30|) valid for convex domains, the signs are 
changed at the odd powers of l/(sR). This means that the result ([317)) is ap- 
plicable also to concave parts of the domain boundary if we take the negative 
value of the outer curvature —1/R(u). 

3.4- Comparison with Savo's results 

It is instructive to compare our coefficients {<7j} of the large-s expansion 
(f30|). expressed via the local curvature in equations (|3T ]) -(f36 l) . with the exact 
ones derived previously by Savo[I3,E|. Savo's results are expressed as integrals 
over the boundary dfl of certain differential operators acting on the shortest 
distance l(r) of a point r from <9f2. For two-dimensional domains, we were able 
to rewrite Savo's complicated formulas in terms of the curvature k(u) and of 
their derivatives. The exact results for the first four coefficients o~i,o~2,o~3 and 
(T4 coincide with our formulas (|3T]) -(|34 |) . In higher orders, we find 

a 5 = ^p/ af2 d.{25^ M -8[ fc V)] 2 }, (44) 
ff 6 = / d"{l3k 5 (uj)+80k(L J )[k'(uj)} 2 +47k 2 (u)k"(uj)}; (45) 

the expressions for the next coefficients are too complicated in Savo's format to 
handle with. Comparing these exact results with the corresponding formulas 
(|35|) and (f36|) we see that the approximation by the local curvature involves 
correctly the powers of the boundary curvature, but misses derivatives of the 
curvature along the boundary. We shall see in the next section that neglecting 
of curvature derivatives has only a small effect on the obtained result. The 
advantage of the approximation by the local curvature (I31l) - (|36[) consists in its 
simplicity, numerical adequacy and availability for extremely high orders. 



4. Pade interpolation 

This section deals with estimates of the low-energy part of the spectrum of 
the (minus) Dirichlet Laplacian. The calculation is based on an interpolation 
between the large-s and small-s expansions of r(s). The convergence of the 
results is good especially for the lowest eigenvalue Ai which is chosen as a test 
of our calculations. 

For large t, the survival probability © can be approximated by a sum of 
first few exponentials with the lowest eigenvalues, S(t) ~ Y^JLi lje~ Xjt - Hence, 
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in the Laplace picture (|12p. we have for general domain 
m 2 

r(s) ~ > 9 J r = d + dis + d 2 s 2 + d 3 s 3 + d±s 4 H , (46) 

e -^i s z + A,- 

3=1 J 

where the even coefficients are given by 



rn 2 



j=i J 



and the odd coefficients vanish d 2 fc+i = 0. The coefficients {d 2 k} are not at 
one's disposal for general domains. For simple domains like the disk or the el- 
lipse, the even coefficients do,d 2 , . . . are, in principle, available to an arbitrary 
order [2^. Within the framework of the m-truncation (|46l) . knowledge of the 
first 2m nonzero coefficients d ,d 2 , ■ ■ ■ ,d 4m _ 2 is sufficient to get very good es- 
timates of the to lowest eigenvalues Xj and the coefficients 7?, by solving the 
corresponding set of 2to nonlinear equations. For example, in the case of the 
unit disk with the exact value Ai = 5.783186 . . ., we obtain Ai = 5.784128 . . . 
for to = 1, Ai = 5.783187 . . . for to = 2, etc. The lowest eigenvalue can be found 
alternatively as lim^oo d 2k /d 2k + 2 = — Ai. 

In the previous section, we derived the large-s expansion of r(s) in the 
approximation by the local curvature: 

3=1 

where the coefficients {&j} are given by (|3"Tj) - (l3"6"l) . Our strategy is to find an 
interpolation scheme between this large-s expansion and the small-s expansion 
(|4"6")l of t(s). The lowest eigenvalue Ai will be deduced directly from this inter- 
polation scheme. 

The interpolation method is based on the following steps: 

(i) Proposal of the Pade interpolation: With respect to the first relation 
in equation (|46|) . the small-s behavior of t(s) can be represented as a rational 
of polynomials P2m-2(s)/Q2m(s) in even powers of s. The large-s expansion of 
t(s) (|48|) contains also odd powers of 1/s. If we try to expand the first relation 
of (|46D in 1/s, the odd powers would be missing. To connect both expansions 
by one interpolation formula, we include also odd powers of s in the rational 
Pim—2(s) /Qim{ s )- The interpolation function is then the Pade approximant 

f\ r / , 01 — Po+PisH hpn-is™" 1 + s™ 

t(s) = \n/n + 2 = n -77, 49) 

W L/ J qo + qis + ■ ■ ■ + q n +is n+1 + s»+ 2 ' V ' 

where for simplification of the notation we use n = 2m — 2 (to obtain more data, 
also odd values of n will be taken). 

(ii) Determination of Pade parameters: To fit the n unknown pj and 
n + 2 unknown qj Pade parameters, we expand the Pade approximant (|49p 
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Table 1. The coefficients of the small- s expansion (46) and the low- 
est imaginary parts of poles of the Pade approximant for the disk 



Pade 


do 


d 2 


di 


rf 6 


Im[s] 


[1/3] 


0.1743 


-0.07472 


0.008383 


-0.00004709 


1.756 


[2/4] 


0.1475 


-0.03538 


0.011763 


-0.00246499 


1.940 


[3/5] 


0.1378 


-0.02803 


0.006355 


-0.00162755 


2.074 


[4/6] 


0.1331 


-0.02516 


0.005090 


-0.00106198 


2.178 


[5/7] 


0.1306 


-0.02371 


0.004541 


-0.00088446 


2.252 


[6/8] 


0.1290 


-0.02287 


0.004246 


-0.00079754 


2.299 


[7/9] 


0.1289 


-0.02235 


0.004067 


-0.00074735 


2.328 



exact 0.1250 -0.02083 0.003581 -0.00061849 2.405 



from both sides s — > oo and s — > 0, in the variables 1/s and s, respectively. 
The expansion in 1/s has to reproduce the series (|48| up to the term of the 
order l/s ra+4 , implying in this way n + 2 conditions for the coefficients pj and 
qj. The expansion in s has to respect the symmetry of the small-s expansion 
(|46p . keeping only even powers of s nonzero. Therefore, we add n conditions 
di = 0, da = 0, • ■ •, G?2n-i = 0. We end up with the complete system of 2n + 2 
nonlinear equations for pj and qj . From all possible solutions we choose the real 
ones. This scheme, which leads to a nonlinear set of equations, differs from the 
standard Pade fitting method, which uses input series data from only one side 
and therefore leads to a linear set of equations. 

(iii) Determination of Ai from Pade approximant: From the first relation 
of equation (|46|) we see that the poles of t(s) occur in the complex plane at the 
points ±i^ < /Xj (j = 1,2,...). The poles of the rational (14*91 are given by 

qo + qis + ■ ■ ■ + q n +i-s n+1 + s n+2 = 0. (50) 

Since the coefficients {qj} are real numbers, the roots of this equation arise in 
complex conjugate pairs. As n — > oo, their (negative) real part should approach 
0. The lowest eigenvalue Ai is determined as the square of the pole closest to 
the origin. 

As a test model, we consider the disk of radius R = 1 for which the exact 
small-s and large-s expansions are given by the formulas (|2l>| and (|2T|) . respec- 
tively. We carried out numerical calculations up to n = 7. In Tab. 1, the first 
few nonzero coefficients of the small-s expansion (|47[) are shown to converge very 
well to their exact values. Fig. 2 depicts the poles of the Pade approximant (|49l) 
which are closest to the imaginary axis and converge to s = ±i\/Ai = ±izo,i for 
n — > oo. It is seen that by increasing the Pade order n, the (always negative) 
real part of the pole goes quickly to zero, as expected. The imaginary part ap- 
proaches the exact value denoted by a cross, with the relative error 0.02 (after 
extrapolation) . The values of these imaginary parts of the poles are presented 
in Tab. 1. In contrast to other numerical methods, like the harmonic inversion 



method 27] or the finite elements method 28], our analytic treatment applies 



easily to general domains. 
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Figure 2: The lowest poles of Pade approximants [n/n + 2] (n = 1, 7) for the disk. 



For the elliptic domain with eccentricity e € (0,1), we complete the in- 
tegrations in (|31j) in polar coordinates, in which the boundary is given by 
r(ip) = bj \J\ — e 2 cos 2 ip, b stands for the minor semiaxis. The dependence 
of the function Ai (e) on the Pade order n is presented in Fig. 3. We see that 
by increasing n our estimates approach the exact plot[29j|. 

The comparison of the results for the function Ai(e), evaluated by using 
the large- s coefficients in the approximation by the curvature (|3ip - (|36[) and the 
exact ones 12. 13| [see (|44| and (|45l) ]. are presented in Fig. 4. The results are 



almost identical up to very large eccentricity e ~ 0.9, the difference remains 
small also beyond this point. 



5. Conclusion 

In this work, we studied the survival probability S(t) of a particle diffusing 
in a finite domain with a smooth absorbing boundary. Its short time expansion 
in © is determined by the boundary characteristics, whilst its long time 
exponential decay is governed by the lowest eigenvalue of the minus Laplacian 
defined in this domain. Thus this quantity enables us to study the impact of 
the boundary shape on the lowest eigenvalue Ai . 

First we proposed a simple method which allows us to calculate the coeffi- 
cients ex,- of the short time expansion of S(t) up to high orders. We approximated 
the boundary locally by circles of the radius inverse to the local curvature of the 
boundary. Using the known exact solution for the disk domain, we expressed 
the coefficients (jj in terms of powers of the curvature integrated along the cir- 
cumference of the domain, see equations (|3ip - (f31)j) . This treatment gives exact 
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results up to 04. In higher orders, the exact coefficients of Savo[12l 1 1 31] contain 
also derivatives of the curvature, which are neglected in the present approxi- 
mation. While Savo's results can be obtained in practice only in the first few 
orders, our formulas are at disposal for virtually any order. 

Next, we showed that there exists an interpolation scheme connecting the 
short time expansion of S(t) with the lowest eigenvalues of the Dirichlet Lapla- 
cian. It is represented by the [n/n + 2] Pade approximation (|49l) of the Laplace 
transform r(s) of the survival probability S(t) (fT2"j) . If the coefficients pi, qj are 
fixed in a specific way, the poles closest to the imaginary axis tend to ±i^fXl. 
Extrapolation of our data for n < 7, calculated for the disk and the elliptic 
domains, exhibits precision within 2% of the exact value. 

Of course, there may exist other interpolation schemes which provide faster 
convergence to the exact value of Ai. Our present aim was mainly to point out 
the existence of such an algorithm. Another task is to extend the calculation of 
S(t) to domains with cusps (with piecewise smooth boundaries), or to higher- 
dimensional domains. These are open problems for future. 
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